Optimal Base Station Placement: A Stochastic Method 
Using Interference Gradient In Downlink Case 
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Abstract 



In this paper, we study the optimal placement and optimal number of base stations added 
to an existing wireless data network through the interference gradient method. This proposed 
method considers a sub-region of the existing wireless data network, hereafter called region 
O , of interest. In this region, the provider wants to increase the network coverage and the users 

throughput. In this aim, the provider needs to determine the optimal number of base stations 
to be added and their optimal placement. 

The proposed approach is based on the Delaunay triangulation of the region of interest and 



I the gradient descent method in each triangle to compute the minimum interference locations. 

^ ' We quantify the increase of coverage and throughput. 

00 

^ 1 Introduction 

• I— I , Our objective is to increase the coverage and the capacity of a network, by adding new base stations. 
^ I Therefore, we need to find the optimal number and optimal placement of additional base stations 
I in an already deployed wireless data network. We consider a sub-region of the existing wireless 

data network, where the service provider wants to increase the coverage of the network and the 

throughput of the users. 

The problem is to find the positions of K additional base stations (a variable number to be deter- 
mined), within the region of interest in an existing wireless data network. Note that the addition 
of K new base stations to the network may impact the coverage and the capacity of the existing 
base stations which makes the problem very difficult. We show that the problem of finding the 
optimal location of a set of new base stations when there is a discrete number of mobile users is 
an NP-Hard problem. The problem in its more general form is much harder when we consider a 
random distribution of the mobile users, even for the uniform distribution of the users. We propose 
a related problem which is to find the points of minimum interference in an existing network. The 
optimality is then to find the minimum interference set of base stations. 
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Figure 1: Atz{/3) of transmitters lying in the central square area of 1000 x 1000 square km. Region 
bounded by a square of side 500 km is the region of interest to place additional base stations. (3 = 1 
and a = 4.0. 

We have studied the optimal placement and optimal number of base stations added to an existing 
wireless data network through the interference gradient method. Our approach is based on the 
Delaunay triangulation of the region of interest (where the triangulation completely covers the 
region of interest). We propose to take as the vertices of this Delaunay triangulation the positions 
of the existing base stations. This is justified by our proof that over each triangle of the considered 
triangulation there is only one global minimum (over the triangle) of the interference function. In 
consequence, there will be as many candidates of minima over the region of interest as triangles in 
the Delaunay triangulation. The problem is reduced to find over this set of candidates the minimum 
interference set of K base stations. We propose two heuristics to find this minimum interference 
points. Numerical simulations show that the coverage and the throughput through this method are 
highly incremented when our method is considered. 

2 Related works 

Plastria pLj presented an overview of the research on locating one or more new facilities in an 
environment where other facilities already exist. The authors of analyze the coverage and the 
capacity of a wireless network regularly distributed. Gabszewicz and Thisse [3] provided another 
general survey on the location problem. Buttazzo and Santambrogio studied the location problem 
but their work does not consider the case when there is an existing wireless data network. Altman 
et al. |S] studied the case when there are two providers in the uplink scenario of a cellular network 
and the users are placed on a line segment. Prommak et al. |6] design a constraint satisfaction 
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problem where they consider the transmitted power, frequency bands, and location constraints. 
The authors of [7j studied a delay-constrained information coverage which targets an optimal base 
station placement with the objective of maximizing information collection within a constraint time, 
which is more suitable for wireless sensor networks. 



We consider an existing wireless data network. In particular, we put ourselves in the context of an 
LTE (OFDMA-based) system, where the relevant objective function is to find the optimal number 
of base stations and their optimal placement to maximize the throughput utility of the data flows. 

We consider OFDMA systems supporting fractional frequency reuse (FFR) for interference mitiga- 
tion similar to the one used on [HI |9] . This type of interference mitigation divide frequency and time 
resources into several resource sets. Fractional frequency reuse in the context of OFDMA systems 
has been discussed in cellular network standarization such as Third Generation Partnership Project 
(3GPP) and Third Generation Partnership Project 2 (3GPP2) [TOl[n]. 

The model is given by the following: 

• Let S = {1, . . . , S} represent the set of base stations which are randomly distributed over the 
two-dimensional plane according to a Poisson point process (PPP) of intensity A. We denote 
by Zj = {xi, yij^ the position of base station i E S. 

• and J sub-bands j E J' = {1, . . . , J} where we denote W the bandwidth of each sub-band. 

• Each sub-band consists of a fixed number c of sub-carriers. 

Furthermore, time is divided into slots consisting of a number of OFDMA symbols and transmissions 
are scheduled to users by assigning a set of sub-carriers on specific slots. The time is slotted, so 
that transmissions within each cell are synchronized, and do not interfere with each other. To 
simplify the exposition we assume that the resource sets span the entire time period. Extension 
to more general resource sets is straightforward [SI [9]. A transmission in a cell, assigned to a sub- 
carrier in a sub-band, causes interference to only those users in other cells that are assigned to the 
same sub-carrier on the corresponding sub-band. We assume that there is no interference between 
sub-carriers. 

• We denote by Tjj(z) G [0,1] the fraction of time an algorithm chooses the user located at 
position z in sub-carrier i in sub-band j. 

• We denote by Rij{z) G [0, B], B < oo the transmission (nominal) rate in sub-carrier i in 
sub-band j, if the user located at position z is chosen. 

Then the average rate a user located at position z actually receives is 



3 The model 




(1) 



We denote the vectors by bold fonts. 
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Each base station may need to allocate a set of sub-carriers and average power P^p[z), if the user 
located at position z is to be assigned to sub-carrier i in sub-band j. 

Let us denote by Gf^\z) the propagation gain from cell to a user located at position z in sub- 
carrier i in sub-band j. 

Then, we consider the SINR (Signal to interference plus noise) function: 

SIN%(z) = ^ ^ , (2) 

where Nq is the thermal noise. 

We use Shannon formula for the rate 

C(z) := W\og^{l + SIN%(z)). (3) 

Let us consider that the propagation gain from cell to a user located at position z associated to 
sub-carrier i in sub-band j is given by the channel gain 



Gf{z):=^=^, (4) 

where dij is the distance between the base station BSj and the mobile terminal MT_,-, h is the height 
of the base station, k is a constant factor, and a > 2 is the path loss exponent. 

To simplify the exposition we consider the distance between base stations to be much greater that 
the height of the base station in eq. (jlj), and k is constant. In the following we will drop the sub- 
indices and analyze just the interference between sub-carriers. This scenario is general and it may 
be applied to any interference problem. 

The objective of this work is to find the optimal positions of K (variable) additional base stations, 
within a bounded region of interest in an existing wireless network, such that it maximizes the 
coverage and the capacity of the existing base stations and these K additional base stations. Note 
that the addition of K new base stations to the network may impact the coverage and the capacity 
of the existing base stations which makes the problem very difficult. 

The coverage of a base station is given by the /3-SINR threshold at which the basic service is 
provided. Considering that P-p is constant for all i, j, k, and assuming negligible thermal noise, we 
define the /3-reception area of the base station i E S, denoted by Aiziif^), which represents the area 
where the transmission from base station i is received with signal to interference plus noise ratio 
(SINR) greater or equal to /3, i.e., 

Amm = {z: |z-zr">/3$^iz-z,r"} (5) 

where |z| represents the Euclidean norm of the vector z = (x, y), i.e., |z| = a/x^ + y'^. 

The coverage of the network is given by the sum over all the base stations of their /3- reception areas 

Coverage(/3) = 5^^7^.(/3). (6) 

ies 



K 
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The /3- reception area of the base station i, AmiP), is computed using the gradient SINR method- 
ology described in [T2] . 



4 Description of our scheme 

4.1 Points of maximum capacity 

Let us assume that the set of S existing base stations is given by 

S = {1,...,S}. (7) 

Let us consider the case we want to incorporate a new base station i ^ S. Then the set of 5 + 1 
base stations will be given by 

5' = 5U{z}, (8) 

The SINR function of a mobile terminal located at position z being served by base station BS^ with 
k E S' located at position z^ is 

SINRfc(z) = = ^, (9) 



where / = |z — Zj|^" and g = — Zj\~°' 

We want to compute the maximum of the capacity function obtained by adding a new base station 
depending on the location of the already existing base stations. The positions of the existing base 
stations are given, and we assume that the distribution of the mobile terminals is uniform, then the 
expected capacity function by adding the new base station is given by 



A 



E(C(z)) = / C{z)dz. (10) 

where 



C(z) = ^ ly log2(l + SINR,(z)). (11) 

The optima of this maximization problem are found at stationary points, where the first derivative 
or the gradient of the objective function is zero. An equation stating that the first derivative equals 
zero at an interior optimum is sometimes called a "first-order condition" . 

We want to study the behavior of the gradient of eq. (ITU]) : 



C(z)rfz . (12) 
\JA J 

We can rewrite eq. (fTTj) as follows: 

^iyiog2(l + SINRfc(z))= (13) 

k£S' 

J2 W^log2(l + SINR,(z)) + iyiog2(l + SINRi(z)) (14) 
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Figure 2: Finding points of local minimum interference in the central square area of 500 x 500 
square km. 



Then in the case when there is a large number of base stations the term 




The relevant term in eq. f lT2|) is then giving by the following 

V / iyiog2(l + SINRi(z))rfz (16) 

J A 

In the low-SINR regime the problem of maximizing the existing capacity is equivalent to the problem 
of maximizing the SINR. In the high-SINR regime both quantities are related however we do not 
show a theoretical result. However, simulation results suggest that even in this scenario we obtain 
a good approximation. In this case, we consider the following function 



A 



V / SINR(z)ciz (17) 

Notice that we can rewrite eq. ([9]) 

SINRfc(^) = ^" = (18) 



|z-z,.|-" 













VSINR(z) = V{^]= J^^ (19) 

\9 J 9 
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In general, this problem is NP-Hard (consider it as the extension of the coverage problem to the 
random distribution of users). Then we propose instead of considering the whole SINR function 
just to consider F{zk) = f{zk,z)dz to be a constant function, then in the previous equation 
we are interested on the behaviour of Vg. 

The problem of maximizing coverage in wireless data networks with K additional base stations 
can be reduced to the classical maximum coverage problem in computer science [13]. In maximum 
coverage problem, the inputs are several given sets and a number K. These sets may have some 
elements in common. The goal is to select at most K of these sets such that the maximum number of 
elements are covered, i.e., the union of the selected sets has maximal size. Formally, the maximum 
coverage problem can be defined as: 

Inputs: A number K and a collection of sets 

S = Si,S2,...,Sm (20) 



Objective: Find a subset 5" C S* of sets, such that |5"| < K and the number of covered elements 
IJ Si is maximized. 

The problem of maximum coverage with given base stations and K additional base stations can 
be reduced to the above described maximum coverage problem. In case of wireless networks, the 
sets Si represent base stations. The elements of each set may represent the users each base station 
may cover with given SINR threshold (3. Note that these sets may overlap for /3 < 1. For /3 > 1, the 
zones covered by the base stations do not overlap [12]. The location of K additional base stations 
shall be chosen so that the coverage of users, in the network, hj N + K base stations is maximum. 
This is analogous to selecting K sets in maximum coverage problem which is an NP-hard problem. 



4.2 Points of minimum interference 

The interference or signal level received at any location z on the two-dimensional plane can be 
represented by the function 

g{z) = Y,\^-z,r (21) 

165 

There can be many local minima and to identify those points, we have used the following two-step 
approach. 

1. We subdivide the region of interest in the network, by using the Delaunay triangulation 
scheme, with the location of the base stations. 

2. We use the gradient descent method to locate the point of minimum interference in each 
triangle. Proof of the convexity of the interference function in each triangle is included in 
Proposition 14.21 
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4.3 Delaunay Triangulation 

A Delaunay triangulation for a set of points S in the two-dimensional plane, denoted by DT(i5) is 
a triangulation (or subdivision of the two-dimensional plane into triangles) such that no point in S 
is inside the circumference of any triangle in DT(»S). 

The Delaunay triangulation of a discrete set of points S corresponds to the dual graph of the Voronoi 
tessellation for S. 

We choose the Delaunay triangulation because we have the following property. 

Proposition 4.1. Any local minimum interference over the triangle T \ {ai, a2, 03} is also a global 
minimum interference over the triangle. 



Proof. Consider the function defined in the two-dimensional plane except in zf. 

g{z,z,):= ^ (22) 
|z - Zi|" 

Then the function of interference is given by 

^(z) = ^^(z,z,) (23) 
ies 

Let us consider the triangle of vertices Oi, 02, 03 where the three vertices are the position of three 
base stations: 

3 3 
r={z : z = ^A,ai;0<Ai<l, {1,2,3}, ^Ai = l} (24) 

1=1 i=l 

Proposition 4.2. The function g{z, A-r,, Zj) is convex in z over T\ {oi, 02, a^}. 

Proof. We just need to prove that ^ is positive. 

Let Tj := |z — Zj|. The gradient of g{z, Zj) in z = {x, y) is: 

d'O I i-a-l Z — Zj Z — Zj , . 

— — —a\z — Zi\ - — —a-, — — . (25 

dz ' ' z-Zj z-Zj"+2 ^ ' 



This means that 



dg x-Xi dg y - yi 

— = —a — To- and — = —a — 
dx r- dy r 



X — xj , y — yi /r.^\ 

— " — -zTo- and — = —a— -To- (26) 



dz2 



= —a 



—a 



1 

|z - Zj|°+2 
1 

|z - Zj|<^+2 



+ (-a-2)|z-Zj 



z — z,- 



- (a + 2)- 



Z - Zj 



Z — Z, 



a{a + 1)- 



z — z, 



|a+2 



(27) 



Since there is no element of S included on T \ {ai, 02, 03}, then the function g is convex on T \ 
{01,02,03} 

The sum of convex functions is convex, then over the domain T \ {fli, 02, 03}, the function 

(?(z) = 5^(7(z,z,) (28) 

is convex over the domain T \ {ai, 02, 03}. 

Any local minimum of a convex function is also a global minimum. 

Important note. Notice that the fact that we are restricting our domain to be inside the triangle 
help us to determine that there is only one global minimum since the domain is a convex subset. 

4.4 Gradient descent method 

Gradient descent method [11] is based on the observation that if the real- valued function is defined 
and differentiable in a neighborhood of a point z°, then the function g{z) decreases fastest if one 
goes from z" in the direction of the negative gradient of g at z°, — V(7(z°). Since our objective is to 
find the minimum of g, then we propose to use the gradient descent method. It follows that, if 

zi - z° - 6t^^^ (29) 
~ l|V^(zO)|| ^2^^ 

and therefore: 

~ l|V^(z")|| 

where 5t > is the step size. Note that g{z'^~^^) < g{z"'). First approximate location of minimum 
interference point, z^, is the centroid of the triangle: 

^0 — yO^ _ ^Pl + ^P2 + ^P2 Upi + yp2 + Vpz ^ ^g-l^^ 

where pi,P2 and are the vertices of the triangle. In order to ensure that the points z" [n > 0) lie 
inside the triangle, we have used the method described in [15] 



4.5 Heuristics 

Here we propose two heuristics to add K additional base stations at the points of minimum inter- 
ference identified by the gradient descent method [T?|. Let V be the set of the location of these 
points of minimum interference. 

1. Heuristic 1: The optimal position of an additional base station is found as follows. 

(a) Rank points in the set V in the ascending order of interference or signal level received 
from all base stations in the set S. 
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Figure 3: 5 base stations, denoted 1 to 5, are added in the region bounded by a square of side 
500 km using Heuristic 1. The number indicates the order of addition to the network. /3 = 1 and 
a = 4.0. 

(b) Select the lowest ranked element p from the set V. The location of p shall become the 
position of the additional base station. 

(c) Remove p from the set V, i.e., V := V — p. 

(d) For additional base stations, repeat from step (b). 

2. Heuristic 2: The addition of a base station at minimum interference point in set V may 
increase the interference at other points of set V located in the immediate vicinity. Therefore, 
remaining points in set V, {V — p}, should be re-ranked. The Delaunay triangulation in 
this case will be provided by the Bowyer- Watson algorithm which gives another approach for 
incremental construction. It gives an alternative to edge flipping for computing the Delaunay 
triangles containing a newly inserted vertex. 

Steps (a) and (b) of Heuristic 2 are the same as in Heuristic 1. After addition of a base station, 
the region of interest is re-triangulated to identify new points of minimum interference. In 
other words, the points of minimum interference in the region of interest shall be identified 
before the addition of all new base stations. 



5 Numerical Simulations 

We consider a wireless data network over an area of 10000 x 10000 square km and base stations 
distributed according to a Poisson point process (PPP) of intensity A = 1 base station per unit 
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Scenario 


Heuristic 1 


Heuristic 2 


Capacity 


1.3559 


1.5614 


1.7007 


Percentage Increase 




15.15% 


25.42% 



Table 1: Average capacity over the surface [bits/s/Hz/km^] over the network and percentage in- 
crease of the capacity under different heuristics 





Scenario 


Heuristic 1 


Heuristic 2 


Total coverage area 


168776 


190501 


204628 


Coverage percentage 


67,51% 


76,20% 


81,85% 


Percentage Increase 




12.87% 


21.25% 



Table 2: Coverage of the network (in square km) and percentage increase of the coverage under 
different heuristics 



square mile. We consider the path loss exponent a = 4. 

In order to ignore edge effects, while computing minimum interference points, we assume that our 
region of interest lies in the center of this network and is a square region with each side of length 
500 km. 

Figure [1] shows the /3-reception areas of the base stations distributed according to a Poisson point 
process of intensity A = 1. It also shows the region of interest marked by a thick lines inside the 
network. 

Figure [2] shows the process of identification of the points of minimum interference in this region 
of interest. Locations marked by 'x', show the first approximation of the points of minimum 
interference. Gradient descent method uses the first approximate position to arrive at the final 
location, marked by 'o', which is the true point of minimum interference in the region formed by 
the Delaunay triangle. Note that, in cases where a triangular region lies partially within the region 
of interest, we are interested in identifying the point of minimum interference lying within the region 
of interest. 

Figures |3] and H] shows the positioning of 5 additional base stations in the network by using the 
approaches of Heuristic 1 and Heuristic 2 respectively. In order to show the improvement achieved 
by the two heuristics, we have computed the total area covered by the base stations with SINR at 
least equal to /3, i.e., Sie^' -^mW) where S' is the set of base stations, also including the additional 
ones in case of Heuristic 1 and Heuristic 2, lying in the region of interest. 

Scenario indicates the existing wireless data network to which we want to increase its coverage 
and its capacity. 

Table 1 summarizes the results obtained by both heuristics. We compute the capacity over the 
network under three different scenarios: Scenario which indicates the capacity on the existing 
wireless data network, and Heuristic 1 and Heuristic 2 explained above. 

Table 2 summarizes the results obtained by both heuristics. We compute the coverage of the network 
under the three different scenarios described above. 
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Figure 4: 5 base stations, shown in red, are added in the region bounded by a square of side 500 km 
using Heuristic 2. The number indicates the order of addition to the network. /3 = 1 and a = 4.0. 

6 Conclusions 

We have studied the optimal placement and optimal number of base stations added to an existing 
wireless data network through the interference gradient method. The problem of finding the optimal 
location of a set of new base stations when there is a discrete number of mobile users is NP-Hard. The 
problem considered is even harder when we consider a random distribution of the mobile users even 
for the uniform distribution. Our proposed method considers a sub-region of the existing wireless 
data network (hereafter called region of interest), where the service provider wants to increase the 
coverage and throughput of the users. Our approach is based on the Delaunay triangulation of the 
region of interest (this triangulation completely covers the region of interest), and it takes as the 
vertices of this triangulation the positions of the existing base stations. We prove that over each 
triangle there is only one global minimum of the interference function. In consequence, there will be 
as many candidates of minima over the region of interest as triangles. Then the problem is reduced 
to find over this set the minimum interference set. We propose two heuristics to find the minimum 
interference points. Numerical simulations show that the coverage and the throughput through this 
method are highly incremented. 

7 Perspectives 

There are many possible extensions to our work, but we think that the more interesting ones are 
the following: 
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• If we want to take into consideration a more realistic LTE (OFDMA) system, the mobile user 
may receive a set of sub-carriers. 

• We have neglected the fact that the distribution of the users may be non-homogeneous over 
the region of interest. 

• In a future work, the shadowing fluctuations over the network will be taken into account. 
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